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ui ; Abstract 
O 

The Schrodinger equation for an electron on the surface of an elliptical torus in the presence of 

a constant azimuthally symmetric magnetic field is developed. The single particle spectrum and 

C/3 . eigenfunctions as a function of magnetic flux through the torus are determined and it is shown 

CO 

^vq that inclusion of the geometric potential is necessary to recover the limiting cases of vertical strip 

and flat ring structures. 
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1. INTRODUCTION 



Toroiclal nanostructures present intriguing possibilites for use as nano-device elements 
In addition to Aharanov-Bohm and persistent azimuthal current effects known 



to exist in quantum rings 



toroidal structures allow for motions around 



the minor radius of the torus subject to boundary conditions distinct from those for flat 
rings (Figs. 1-3). 

For hollow torii, electrons are thought to be localized near the surface of the object. The 



restriction to motion near a surface has interesting manifestations. Recent work 
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on 



a torus of major radius R and circular cross section of radius a, has shown that a surface 



dependent geometric potential Vc 



Q 15, 16, Q 18, Q 



is important even for electrons 



that can wander substantial distances away from the surface and should be employed as an 
effective potential when considering two dimensional problems on curved surfaces. Curvature 
effects are more pronounced for elliptical torii (ET 2 ) which can behave in different ways than 
a torus with circular cross section. A recent paper by Gravesen, Willatzen and Lew Yan 
Voon in this journal |20| dealt with particles constrained to motion on surfaces of revolution, 
including elliptical torii. There the influence of toroidal eccentricity and Vc were shown to 
strongly affect system eigenvalues and eigenfunctions. This work is concerned with the 
extension of |20j to include a constant magnetic field along the z-axis, and the role Vc plays 
in recovering the ribbon (taken to be an infinitely thin vertical strip of radius R and height 
2b) and flat ring (an annular region of inner radius R — a and outer radius R + a) limits. 

This paper is organized as follows: section 2 presents the geometry leading to, and 
the formalism by which, the Hamiltonian on ET 2 is derived inclusive of a magnetic field 
B = Boe z . The methodology presented in section 2 differs from that given in 20] because 
here the magnetic field is to be incorporated via the minimal prescription so it proves 
advantageous to first derive the gradient operator rather than directly employ 



V 2 — ~^ ^ 



9 2 9 



_d_ 



(1) 



In section 3 numerical results are given as curves of single particle ground state energies 
Eq as functions of magnetic flux 7 for elliptical torii of several eccentricities. Section 4 is 
reserved for conclusions and suggestions for future work. 
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2. DEVELOPMENT OF THE HAMILTON! AN ON ET 2 



As noted above, including a vector potential A in the Hamiltonian makes it advantageous 
to derive the gradient operator and proceed to include A by the minimal prescription. 

Let (e p , e<£, e 2 ) be cylindrical coordinate system unit vectors. Parameterize points near 
an elliptical toroidal surface of major radius R and minor radii a and b by 

x(0, 0, q) = (R + a cos8)e p + b sin#e 2 + qe n (2) 

with e n the unit normal to ET 2 to be defined momentarily and q the coordinate measuring 
the distance from the surface. In what follows the major radius R will be set to R = 500A, 
which is in accordance with fabricated structures jj, 21, 2^|. From Eq. (2) 

dx. = PdOei + Wd(f)e^ + dqe n + qde n (3) 

with 

P = (a 2 sm 2 9 + b 2 cos 2 6) 1/2 (4) 



W = 1 + a cos# (5) 



e x = —{—a sin#e p + b cos#e 2 ) (6) 

and 

e n = — (a sin^e^ + b cos9e p ). (7) 

The unit vectors (ei,e n ) are tangent to ET 2 in the direction of increasing 9 and <p respec- 
tively. 
From 

de n = eid9 + cosOe^dcj), (8) 

the vector differential line element can be found from Eqs. (8) and (3), leading to a gradient 
operator 

- 1 d_ 1 d_ d_ 

~ P(l + k e q) ei m + Wil + k^qf^ + [ ) 

with kg, k<j, the principle curvatures given by 

, ab b cos9 . . 

fcfl=p3, k <t> = -Wp- (10) 



The vector potential (working in the Coulomb V • A = gauge) appropriate to B e z is 

A=f(nr + ?^!*)„. (11) 

The Schrodinger equation (h = m = 1) 

i^lv + gA^ \P = (12) 

which results from Eqs. (9), (11) and (12) can be reduced to a concise and dimensionless 
form by first 



a. noting 



<9ei 

~de 



de r 



ab 
~p2 en 



ab 



09 P 2 



dei a sin9 



b cosd 



(13) 
(14) 
(15) 
(16) 



b. setting a = a/R, j3 = b/R, e = 2Ea 2 , 

c. letting D(0) = P(0)/R, p(0) = P(6)/a and F{9) = l + a cosO, 

d. defining 7 = .2Q3B , with B in Teslas, which is the conversion factor for an 
R = 500A torus, 

e. performing the well-known procedure for obtaining Vc (which will appear below as 
the scaled dimensionless function Up ) for which the reader is directed to the relevant 



references LL4L LL5 
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28, 



29 



30]. 



f. noting that the azimuthal symmetry of the problem allows for the eigenfunction on 
ET 2 to be taken as ^(6, <fi) = ip(8)exp [iv<fi\. 

Applying the conventions and procedures listed above results in a pair of equations for the 
surface and normal variables in the q —> limit 



d 2 ijj ( a sin 9 a 2 — (3 2 \ dtp 
W 2 ~\ F{9) + D 2 (9) J~d9 



D 2 (9)v 2 



_ F 2 {9) 



+ U C {9)1 
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2 F 2 (9)a 2 p 2 (9) 



^- r yh'a 2 p( y 9) 



V = o, 

(17) 
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-\^ + V n {q) X n = E nXn . (18) 
The scaled dimensionless curvature potential Uq appearing in Eq. (17) is 

um = JJW + ?~*\ (19) 



A\D A (9) F 2 (9) 

The normal confining potential V n (q) can be chosen to take any convenient form, and it is 
apparent in the q — > limit that the surface Schrodinger equation is independent of the 
choice. The independence of surface observables on V n (q) has also been shown to be a good 
approximation when the particle is allowed to move in a finite thin layer by a basis set 
calculation in the full three dimensional space The results that follow will not include 

n 

E n ; it was also shown in |13J there is negligible state mixing in the q degree of freedom so 
that the surface spectrum is essentially independent of E n . 



3. NUMERICAL RESULTS 



The Schrodinger equation given by Eq. (17) is invariant under 9 — > —9 making it pos- 
sible to separate its eigenfunctions into even and odd ^-parity states by proper choice of 
initial conditions at 9 = jsi] ]. Numerically it proves convenient to generate two linearly 
independent solutions of Eq. (17) with initial conditions 

V>(0) = 1, V'(0) = -»• i> A (20) 

V(0) = 0, V ; (0) = 1 -»• ipB (21) 

and insisting upon 

Aip A (0) + Bij B (0) = A^ a {2ti) + Bi{> B (2ir) (22) 



A^' A (0) + Bf B (0) = Aj>' A (2n) + B^(2tt). (23) 

Eqs. (22) and (23) can be rearranged into a homogeneous linear system for the coefficients 
A, B and the entire spectrum quickly determined by a simple loop over e appearing in Eq. 
(17). 
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Single particle ground state energy e plots as a function of 7 are given for three toroidal 
eccentricities in Figs. (4-6). The most interesting feature emerging from inspection of Figs. 
(4) and (5) is the necessity of including Uq to approach the ring and ribbon limit. However, 
as shown in Figs. (7) and (8), the ring and ribbon limits do differ in detail when compared 
to the curves for the elliptical torii; because persistent currents are sensitive to the shape 
of the free energy as a function of flux, the differences in detail may prove important upon 
a more extensive calculation of spectra. Comparison of the curves shown in Fig. (6) with 
the values of the limiting cases appearing in the previous two figures indicate clearly that a 
torus with circular cross-section cannot realistically be approximated with a ring or ribbon, 
and that the Uc alters the structure of £0(7) substantially. 



4. CONCLUSIONS 



In this work single-particle ground state energies Eo as a function of magnetic flux 7 
were calculated for three elliptical torii of very different character. The curves indicate 
that any attempt to model toroidal structures by approximating them as two-dimensional 
ribbons or rings must include the geometric potential in order to recover the respective 
two dimensional limits. This result is surprising in that Vc is large only near the regions 
of substantial curvature and is negligible over most of the structure. Nevertheless, it has 
considerable effect on ET 2 eigenvalues, and its omission leads to disagreement with the flat 
limits. The inclusion of Vc, however, does not simply trivially reproduce the flat limits for 
the toroidal eccentricities investigated here; the £0(7) curves show differing peak heights and 
locations than those of the limiting cases. 

An interesting natural extension of this work is the addition of an off-axis component of 
the applied magnetic field. Those cases are perhaps best suited to a basis set method as 
that employed in subject to modification of the integration measure arising from toroidal 
eccentricity. The interplay of the magnetic field with regions of substantial curvature may 
yield interesting mixing of azimuthal modes depending on the magnetic field orientation. 



6 



5. ACKNOWLEDGEMENTS 

M.E. would like to thank B. Etemadi for useful discussions. 



[1] H.R.Shea, R. Martel, and P. Avouris, Phys. Rev. Lett. 84, 4441 (2000). 

[2] K. Sasaki, Y. Kawazoe, and R. Saito, Phys. Rev. A 321, 369 (2004). 

[3] M. Sano, A. Kamino, J. Okamura, and S. Shinkai, Science 293, 1299 (2003). 

[4] H. Zhang, S. W. Chung, and C. A. Mirkin, Nano. Lett. 3, 43 (2003). 

[5] I. Filikhin, E. Deyneka, and B. Vlahovic, Modelling Simul. Mater. Sci. Eng. 12, 1121 (2004). 

[6] S. Gylfadottir, M. Nita, V. Gudmundsson, and A. Manolescu, Phys. E 27, 209 (2005). 

[7] S. Viefers, P. Koskinen, P. S. Deo, and M. Manninen, cond-mat /0310064. 

[8] S. Latil, S. Roche, and A. Rubio, Phys. Rev. B 67, 165420 (2003). 

[9] K. Sasaki and Y. Kawazoe, Prog. Theo. Phys. 112, 369 (2004). 

[10] D. Gridin and A. Adamou, Phys. Rev. B, volume = 69, pages = "155317" (2004). 

[11] Y. V. Pershin and C. Piermarocchi, cond-mat/0506187. 

[12] M. Encinosa and L.Mott, Phys. Rev. A 68, 014102 (2003). 

[13] M. Encinosa, L.Mott, and B. Etemadi, Phys. Scr. 72, 13 (2005). 

[14] H. Jensen and H. Koppe, Ann. of Phys. 63, 586 (1971). 

[15] R. C. T. da Costa, Phys. Rev. A 23, 1982 (1981). 

[16] R. C. T. da Costa, Phys. Rev. A 25, 2893 (1982). 

[17] P. Exner and P. Seba, J. Math. Phys. 30, 2574 (1989). 

[18] S. Matusani, J. Phys. Soc. Jap. 61, 55 (1991). 

[19] M. Burgess and B. Jensen, Phys. Rev. A 48, 1861 (1993). 

[20] J. Gravesen, M. Willatzen, and L. L. Y. Voon, Phys. Scr. 72, 105 (2005). 

[21] A. Lorke, R. J. Luyken, A. O. Govorov, and J. P. Kotthaus, Phys. Rev. Lett 84, 2223 (2000). 

[22] J. M. Garca, G. Medeiros-Ribeiro, K. Schmidt, T. Ngo, J. L. Feng, A. Lorke, J. Kotthaus, 

and P. M. Petroff, App. Phys. Lett. 71, 2014 (1997). 

[23] P. Duclos and P. Exner, Rev. Math. Phys. 7, 73 (1995). 

[24] J. Londergan, J. Carini, and D. Murdock, Binding and scattering in two dimensional systems; 
applications to quantum wires, waveguides, and photonic crystals ( Springer- Verlag, Berlin, 



7 



1999). 

[25] J. Goldstone and R. L. Jaffe, Phys. Rev. B 45, 14100 (1991). 

[26] P. Ouyang, V. Mohta, and R. L. Jaffe, Ann. of Phys. 275, 297 (1998). 

[27] I. Popov, Phys. Lett. A 269, 148 (2000). 

[28] S. Midgley and J. Wang, Aus. J. Phys. 53, 77 (2000). 

[29] I. J. Clark and A. J. Bracken, J. Phys. A 29, 4527 (1996). 

[30] P. C. Schuster and R. L. Jaffe, Ann. Phys. 307, 132 (2003). 

[31] M. Encinosa and B. Etemadi, Found. Phys. Lett. 16, 403 (2003). 

[32] M. Encinosa, Phys. E 28, 209 (2005). 



8 



Figure Captions 



Figs. 1-3: Bohmian trajectories corresponding to a two eigenfunctions superposition on an 
R = 500A, a = 250A torus with a circular cross section as a function of magnetic field values 
(top to bottom) B = 0, 2, 4 Tesla. 

Fig. 4: 5 for an a — .5,(3 — .1 elliptical torus plotted as a function of magnetic flux 
7 = .263£> - Diamonds correspond to the Uc — case, stars to Uq ^ 0, and squares to a 
two dimensional annular region with inner radius 1 — a and outer radius 1 + a. 

Fig. 5: Sq for an a = = .5 elliptical torus plotted as a function of magnetic flux 
7 = .263£>o- Diamonds correspond to the Uc = case, stars to Uc ^ 0, and squares to a 
two dimensional vertical strip/ribbon with unit radius and height 2f3. 

Fig. 6: e for an a = .5 circular torus plotted as a function of magnetic flux 7 = .263fi - 
Diamonds correspond to the Uc — case, stars to £/c 7^ 0. 

Fig. 7: Detailed plot of e for an a — .5, f3 — .1 elliptical torus plotted as a function of 
magnetic flux 7 = .263i?o- Stars correspond to the limit of the two-dimensional annular 
region with inner and outer radii 1 =F a and diamonds to the elliptical torus with C/ c ^0. 

Fig. 8: Detailed plot of e for an a — .1,(3 — .5 elliptical torus plotted as a function of 
magnetic flux 7 = .263£>o- Stars correspond to the limit of the two-dimensional vertical 
strip/ribbon of unit radius and height 2(3, and diamonds to the elliptical torus with Uc ^ 0. 
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